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TECHNICAL  NOTE  NO.  1875 


ON  TWO-DIMENSIONAL  SUPERSONIC  FLOWS 
Fjy  Stefan  Bergman 


SUMMARY 


In  previous  reports  "by  the  author,  the  method  of  operators  has  been 
applied  to  the  investigation  of  two-dimensional  irrotational  flow  patterns 
of  an  ideal  compressible  fluid,  mostly  in  the  subsonic  case.  Using  a 
similar  approach  in  the  present  paper,  the  author  defines  two  different 
operators  which  generate  supersonic  two-dimensional  flow  patterns  from 
differentiable  functions  of  one  real  variable.  (An  operator  is  any  rule 
by  means  of  which  one  function  is  converted  into  another.  This  concept 
is  discussed  in  some  detail  in  the  body  of  this  paper.)  Since  operators 
of  the  type  considered  preserve  many  properties  of  the  functions  to  which 
they  are  applied  and  since  the  theory  of  functions  of  one  real  variable 
is  more  extensively  developed  than  that  of  solutions  of  the  compressi¬ 
bility  equations,  the  results  obtained  can  be  used  as  a  basis  for  the 
investigation  of  supersonic  flow  patterns. 


INTRODUCTION 


Ify  following  the  line  of  approach  developed  in  previous  publications 
for  the  study  of  two-dimensional  subsonic  flows  (see  references  1  to  4), 
the  method  of  operators  is  applied  in  the  present  paper  for  generating 
two-dimensional  steady  supersonic  flow  patterns.  Using  the  holograph 
method  introduced  by  Chaplygin,  the  author  considers  the  stream  function 
in  the  A, 9— plane,  where 


A  =  h  ■'•arc  tan  jh(M^  —  l)-'-/^  —  arc  tan  (vfi  —  l)-'-/^J 
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(see  equations  (5)  and  (6)) ,  M  teing  the  local  Mach  number;  Q}  the 
angle  "between  the  velocity  vector  and  the  x— axis;  and  h,  a  constant 
characteristic  of  the  fluid.  ^ 

The  stream  functions,  considered  as  a  function  of  A  and  0 , 
satisfy  a  linear  differential  equation 


*AA  -  *00  +  4W1*A  =  0 


Two  operators  which  generate  solutions  of  this  equation  are  defined 
in  the  present  paper.  The  first  operator  is  obtained  by  using  Kiemann's 
function.  (See  section  entitled  "Operator  Obtained  by  the  Use  of 
Kiemannfs  Function. " )  The  second  operator  is  derived  by  means  of  the 
theory  of  integral  operators.  In  the  section  dealing  with  the  second 

operator,  a  set  of  functions  H(A),  E^fA),  n  =  1,  2,  .  .  . ,  are 
determined  such  that 


00 

] 

P(f)  =  lim  - 

H(A) 

f(D  +  y~  Em(n)(A)ftn](o 

m — » 00 

n=l 

J 

^At  every  point  (x,y)  of  a  steady  two-dimensional  flow  the 
velocity  vector,  the  Cartesian  components  of  which  are  q  cos  0 
and  q  sin  0,  is  defined.  If  i(x,y)  denotes  the  stream  function 
and  p(x,y),  the  density,  the  coordinates  x  and  y  can  be  expressed 
as  functions  of  q  and  0  (and  hence,  using  relations  (5)  and  (6) 
of  A  and  0)  by  solving  the  equations 


q  cos  0 


-1  dt 
P  *r 


and 


q  sin  0 


Substituting  these  expressions  for  x  and  y  into  t(x,y),  the 
function  Y(a,0)  =  t[x(A,0),  y(A,0)]  is  obtained,  i.e.,  the  stream  function 
in  the  A, 0— plane. 
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where 


f^(£)  =  f^n  ^(1)  d! 


f[°](l)  =  f(6) 


£  =  A  +  0  or  |  =  A  —  9 


(where  f  is  an  arbitrary  differentiable  function  of  one  real  variable  | ) 
represents  a  solution  of  the  equation 


*AA  ~  *00  +  4N1*A  = 


so  that  P(f)  can  be  interpreted  as  the  stream  function  of  a  supersonic 
flow.  The  equation  of  a  streamline  of  the  flow  in  the  physical  plane  is 
then  given  in  parametric  form: 


1  .  dA  1  ,  .  _ 

—  tA  cos  0  - - Bin  0 )  ae 

P  A  dq  pq  y  ) 


+  - - - —  to  cos  0  — - —  tA  sin  e  dA 

no2  6  dA  pq  A 


if -  tA  sin  0  M  +  1_  j,.  cos  0  j  d0 
J\P  A  dq  pq  ^0  } 


(m2  -  l)  dq  1 

- -  \|rfl  sin  0  —  +  —  \|r.  cos  0  dA 


2  e  dA  pqYA 


where  ^(A,#)  =  c  =  Constant.  The  representation  P(f)  is  investigated 
in  the  present  paper ,  and,  in  particular,  methods  for  the  actual 
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determination  of  E^11^  are  discussed  as  well  as  the  domain  in  the 
A ,,0— plane  where  the  order  of  summation  and  passage  to  the  limit  can  "be 
inverted . 

The  evaluation  of  the  operators,  that  is,  the  determination  of  the 
flow  pattern  when  the  function  f  is  given,  involves  lengthy  computa¬ 
tions.  However,  in  a  manner  similar  to  that  described  in  reference  4 
for  the  subsonic  case,  it  is  possible  to  prepare,  once  and  for  all, 
tables  which  greatly  facilitate  numerical  computation  of  flow  patterns. 

The  results  obtained  in  the  present  paper  can  be  used  as  a  basis 
for  the  investigation  of  supersonic  flows. 

The  next  problem  which  arises  in  this  connection  consists  in  deter¬ 
mining  for  a  desired  flow  the  functions  of  one  variable  which  have  to  be 
inserted  into  the  operator  in  order  to  obtain  the  desired  boundary  shape 
or  to  decide  that  such  a  flow  does  not  exist.  The  author  expects  to 
treat  these  problems  in  subsequent  publications.  See  also  the  SUMMARY 
REMARKS  presented  at  the  end  of  the  present  paper. 

This  investigation  was  conducted  at  Brown  University  under  the 
sponsorship  and  with  the  financial  assistance  of  the  National  Advisory 
Committee  for  Aeronautics. 

The  author  would  like  to  express  his  sincere  appreciation  for  the 
assistance  and  advice  he  received  from  Mr.  Bernard  Epstein  and  to  thank 
Mr.  Maurice  H.  Slud  for  his  interest  and  valuable  help  in  connection  with 
the  present  investigation. 

SYMBOLS 


a 


k 

<1 

u 

v 


velocity  of  sound  in  compressible  fluid 
velocity  of  sound  at  a  stagnation  point 


ratio  of  specific  heat  at  constant  pressure  to  specific 
heat  at  constant  volume  (k  =  1.4  for  air) 


speed  of  a  fluid  particle  VVu2  +  v^ j 
x— component  of  velocity  vector  (q  cos  9 ) 
y— component  of  velocity  vector  (q  sin  0) 


Cartesian  coordinates  in  physical  plane 
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B  =  VM^  -  1  for  M  >  1 


M  Mach  number  (q/a) 

R(!  jT);I0Jti0)  Biemann’s  function  for  equation  (15) 


T  =  Nil  —  ¥r  for  M  <  1 


3  =  T  +  iB 


l  =  X  +  iA 


(t) 

Fn  (2 A) 


functions  defined  hy  recurrence  relations  (22) 
see  equation  (l6) 


F2(2A' ) 


see  Theorem  under  section  Integral  Operator  of  Type  Given 
By  Formula  (ll) 


Fg^^A')  a  sequence  of  functions  converging  to  F2(2A’)j  see  equa¬ 
tion  (ko) 

coefficients  of  Taylor  series  expansion  of  Fg(2A* );  see 
equation  (39) 

£  =  A  +  0 
T)  =  A  —  & 

I*.  =  a*  +  e 

V  =  A*  -  9 

X  =  -  (log  +  -  log  - - —  )  for  M  <  1;  log  =  loge 

2  y  .  !  +  T  h  h-1  -  T / 


A  =  ^  arc  tan  (hB)  —  arc  tan  B  for  M  >  1 


A»  =  A  —  —  jt  (  —  —  1 


2  Vh 


6  angle  between  velocity  vector  and  positive  direction  of 

x— axis  (arc  tan  v/u) 


potential  function 
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a* 

* 

stream  function 

* 

>|f* 

"reduced  stream  function"  (see  equation  (l4)) 

a  sequence  of  functions  converging  to  see  equation  (37) 

po 

density  ’of  fluid  at  a  stagnation  point 

rv 

P 

density  of  fluid  in  motion 

P 

"reduced  density"  ^p/pQ^ 

% 

see  equation  (9) 

Note.-  Frequently  quantities  are  considered  as  functions  of  different 
pairs  of  variables,  so  that  different  symbols  should  he  used  to  designate 
the  functional  dependence  in  the  various  planes.  However,  this  is  not 
done  herein,  the  same  symbol  being  used  in  each  plane.  This  should  cause 
no  confusion,  for  the  meaning  is  clear  from  the  context. 

r 

ANALYSIS 

The  Mathemati  cal  Problem  in  Simplified  Form  in  the  Case 

r 

of  a  Supersonic  Motion  of  a  Compressible  Fluid 

In  this  report,  operator  methods  are  applied  to  the  problem  of 
constructing  two-dimensional,  irrotational,  steady  flow  patterns  of  a 
compressible  fluid  at  supersonic  velocities. 

Just  as  in  the  case  of  an  incompressible  fluid,  the  steady,  irrota— 
tional  flow  is  completely  described  by  either  the  velocity  potential  0(x,y ) 
or  the  stream  function  t(x,y) .  If  it  is  assumed  that  the  motion  is  adia¬ 
batic  and  that  the  fluid  is  a  polytropic  gas,  then  the  motion  is  governed 
by  the  following  system  of  nonlinear  partial  differential  equations: 

< 

II 

il 

II 

1 

>  (1) 

f 

8x 
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and 


1 


Here  p  =  p /pQ  is  the  "reduced  density"  of  the  fluids  that  is,  the 
density  p  divided  h y  the  density  pQ  at  rest;  x,y  are  Cartesian 

coordinates  in  the  physical  plane,  that  is,  in  the  plane  where  the 
motion  takes  place;  k  is  the  ratio  of  the  specific  heat  at  constant 
pressure  to  the  specific  heat  at  constant  volume;  and  a0  is  the 
velocity  of  sound  in  the  fluid  at  rest.  (See  reference  1,  equations  (24) 
and  (25).) 

A  fluid  motion  satisfying  equation  (l)  may  show  two  essentially 
distinct  types  of  behavior,  according  as  the  speed  of  the  fluid  is 
greater  or  less  than  the  speed  of  sound  in  the  fluid.  '  The  mathematical 
reason  is  readily  seen  in  the  following  manner.  If  p  and  '  t  are 
eliminated  from  equation  (l)  the  result  is  a  single,  nonlinear,  partial 
differential  equation  of  the  second  order  for  0: 

(a2  -  0/)^  -  +  (a2  -  0y2>„  -  0  (2) 


where 


a2  =  ao  2-|(k-1)fc2  +  0y2)  (3) 

is  the  square  of  the  local  velocity  of  sound.  (See  equation  (28)  of 
reference  1.)  The  discriminant  of  equation  (2)  is 

-  (a2  -  ^%2  -  0y2) 

=  -a2(a2  -  0X2  -  0y2)  =  a2(q2  -  a2)  (4) 

where  q^  is  the  speed  at  the  point  x,y  under  consideration.  In  regions 
where  the  velocity  is  subsonic,  that  is,  q  <  a,  this  equation  is  of 
elliptic  type,  whereas  in  regions  where  the  velocity  is  supersonic,  that 
is,  q  >  a,  equation  (2)  is  of  hyperbolic  type.  (See  reference  5*  pp.  1— 4.) 
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The  solutions  of  these  two  types  of  differential  equation  have  markedly 
different  functional  character,  and  these  differences  are,  of  course, 
reflected  as  the  two  different  types  of  compressible  fluid  motion. 

In  spite  of  these  differences,  however,  the  operators  which  have 
been  used  for  the  solution  of  subsonic  flow  problems  may  be  adapted  to 
the  solution  of  supersonic  problems  after  certain  formal  modifications. 
Before  this  is  done,  however,  it  is  necessary  to  make  a  few  preliminary 
remarks . 

The  differential  equations  (l),  both  for  subsonic  and  supersonic 
velocities,  were  derived  under  rather  restrictive  assumptions,  namely, 
that  the  flow  is  two-dimensional,  that  the  fluid  is  inviscid,  and  that 
the  motion  is  irrotational  and  steady.  By  these  assumptions  the  problem 
is  considerably  idealized;  but  experience  has  shown  that  in  a  great  many 
important  cases  occurring  in  practice,  when  the  speed  is  small,  the 
consequences  of  the  theory  are  in  satisfactory  agreement  with  observa¬ 
tions.  Moreover,  a  physical  justification  for  the  apparently  drastic 
assumption  that  the  motion  is  laminar  is  contained  in  the  Von  Mises 
"hydraulic  hypothesis."  (See  reference  6,  p.  81.) 

In  many  instances  where  the  indicated  idealization  is  not  admissible, 
the  "idealized"  theory  can  at  least  be  used  as  a  first  approximation, 
which  can  be  improved  by  correction  terms.  So,  for  instance,  the 
boundary-layer  theory  describes  the  motion  near  the  boundary  where 
viscosity  effects  cannot  be  neglected,  whereas  in  the  interior  of  the 
flow  pattern  the  theory  of  ideal,  inviscid  fluid  can  be  successfully 
used  without  any  alteration  in  the  case  of  a  comparatively  slow  motion. 

Unfortunately,  when  the  speed  increases  considerably,  the  situation 
changes  completely,  and  the  theory  of  an  incompressible  fluid  can  no 
longer  be  used  as  a  first  approximation.  Furthermore,  at  supersonic 
velocities  a  completely  new  phenomenon  appears,  the  influence  of  which 
is  not  easily  assessed.  When  a  compressible  fluid  flows  around  an 
immersed  body  at  supersonic  velocities,  a  very  characteristic  phenomenon 
is  the  appearance  of  a  shock  line.  It  appears  as  a  curve  (or  curves) 
which  divides  the  physical  plane  into  two' (or  more)  unconnected  infinite 
regions.  In  the  neighborhood  of  the  shock  line,  the  fluid  may  no  longer 
in  general  be  regarded  as  nonviscous;  there  is  a  discontinuity  of 
velocity  along  the  shock  line;  and  in  the  region  of  the  physical  plane 
beyond  the  shock  line  (except  in  a  few  special  cases)  the  flow  is  no 
longer  irrotational.  For  all  these  reasons,  the  differential  equa¬ 
tions  (l)  may  no  longer  be  supposed  to  represent  the  actual  physical 
situation.  Exact  equations  can  be  derived,  of  course,  but  they  are 
extremely  complicated  nonlinear  partial  differential  equations  for  which 
no  mathematical  theory  now  exists;  and  the  problem  presented  (which  is 
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neither  an  initial— value  problem  nor  a  boundary— value  problem,  because 
the  shock  line  is  not  prescribed)  cannot  be  handled  by  any  known  mathe — 
mat i cal  means. 


Shock 

line 


An  example  of  a  supersonic  flow  pattern 
around  an  obstacle  with  a  shock  line. 


It  can  be  expected  that  in  the  same  manner  as  the  idealized  theory 
of  an  incompressible  fluid  represents  an  excellent  "first  approximation" 
in  the  sense  just  explained  for  comparatively  slow  motions,  the  idealized 
theory  of  compressible  fluid  will  represent  a  similar  "first  approximation" 

for  motions  with  considerable  speed2  (if  necessary,  with  certain  addi¬ 
tional  alterations). 

In  this  and  previous  papers  on  the  subject  by  the  author,  therefore, 
no  attempt  is  made  to  find  exact  solutions  of  the  actual  physical  problem. 
Instead  it  is  assumed  that  the  part  of  the  flow  around  the  obstacle 
satisfies  the  conditions  under  which  equations  (l)  were  derived.  The 
method  of  operators  is  then  used  to  find  and  investigate  exact  solutions 
of  this  idealized  problem.  The  results  may  be  regarded  as  first  approxi¬ 
mations  to  the  true  facts  in  the  sense  just  explained. 

2 

It  is  clear  that  for  many  purposes  it  will  be  necessary  to  develop 
a  three-dimensional  theory.  This  theory,  even  in  the  case  of  an  incom¬ 
pressible  fluid,  is  only  in  the  very  first  stages  of  development.  The 
theory  of  operators  and  some  other  modern  mathematical  means,  as  the 
author  expects  to  show  in  another  paper,  will  make  it  possible  to  deal 
with  certain  problems  in  the  three-dimensional  case.  These  methods, 
besides  yielding  results  for  subsonic  motions  of  a  compressible  fluid, 
will  also  yield  results  of  certain  problems  (as  yet  unsolved)  for  an  *■ 
incompressible  fluid. 
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The  question  may  also  he  regarded  from  another  viewpoint.  It  would 
often  he  possible  to  find  a  numerical  solution  of  a  specific  problem,  hut 
such  a  result  has  little  theoretical  value.  A  satisfactory  theory  must 
he  capable  of  doing  much  more.  It  must  he  able  to  give  at  least  qualita¬ 
tive  information  concerning  how  the  motion  of  the  fluid  is  affected  by 
variations  in  the  form  of  the  immersed  body,  the  specific  properties  of 
the  fluid,  or  any  other  significant  variables;  and  preferably  it  should 
also  be  capable  of  providing  quantitative  predictions.  The  operator 
method  of  solving  differential  equations  has  the  advantage  that  it  may 
be  used  to  investigate  the  functional  properties  of  the  solutions  it 
yields.  (For  details,  see  references  2  and  3*) 

Remark.—  From  the  purely  mathematical  point  of  view,  an  operator  is 
a  rule  which  is  applied  to  a  function  of  an  appropriate  class  (usually 
denoted  as  class  A)  to  obtain  a  function  belonging  to  another  class 
(class  C).  A  simple  example  of  an  operator  is  the  process  of  obtaining 
a  harmonic  function  of  two  real  variables  (class  C)  by  taking  the  real 
(or  imaginary)  part  of  an  analytic  function  of  a  complex  variable 
(class  A).  This  operator  preserves  many  properties  of  analytic  functions 
and  is  often  used  as  a  tool  for  investigation  of  harmonic  functions. 

A  similar  relationship  may  be  seen  in  the  operators  employed  in  the 
theory  of  compressible  fluid.  These  operators  appear  in  the  form  (equa¬ 
tion  (11))  for  subsonic  flows  and  in  a  similar  form  for  supersonic  flows. 
In  the  case  of  subsonic  flows,  equation  (ll)  acts  on  analytic  functions 
of  a  complex  variable  (class  A)  and  produces  complex  solutions  of  equa¬ 
tion  (10)  (class  C).  In  the  supersonic  case,  the  operator  acts  on  a 
differentiable  function  of  a  real  variable  (class  A’)  and  produces  solu¬ 
tions  of  equation  (9)  (class  C').  These  operators  are  special  cases  of 
operators  introduced  by  the  author  which  convert  analytic  functions  of  a 
complex  variable  into  complex  solutions  of  a  given  linear  differential 
equation  of  elliptic  type  and  differentiable  functions  of  a  real  variable 
into  solutions  of  a  given  equation  of  hyperbolic  type. 


The  theory  of  integral  operators  gives  a  deep  insight  into  the  property 
of  the  solutions  of  the  compressibility  equation.  As  indicated  in 
references  1,  2,  and  3,  solutions  of  equation  (10)  can  be  written  in  the 
form 


(See  equation  (55)  of  reference  1.)  Here  \(M)  is  a  function  of  the 
Mach  number  M,  which  in  the  subsonic  case  is  real,  and  in  the  supersonic 

case  purely  imaginary.  The  expression  E(M,0,t)  =  E(A.,0,t)  is  a  function 
which,  in  the  "simplified  case"  (Tricomi  case),  when  considered  as  a 

function  of  u  =  t2(X  +  i0)/2A,,  satisfies  the  hypergeometric  equation. 

In  the  case  of  the  "exact"  compressibility  equation,  behavior  and  proper— 

r\j 

ties  of  E(A,,0jt)  are  completely  analogous  to  those  in  the  simplified 
case.  It  is  described  in  sections  k  and  5  of  reference  3. 
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The  operator  method^  has  another  advantage  vhich  is  very  important 
from  the  practical  point  of  view.  All  too  frequently  the  only  thing 
which  stands  in  the  way  of  solving  an  engineering  problem  mathematically 
is  the  difficulty  of  performing  the  incidental  numerical  computations. 
Although  the  use  of  modern  automatic  computing  machinery  in  many  instances 
reduces  considerably  the  amount  of  labor  needed,  it  still  remains  a 
serious  obstacle.  Whenever  possible,  the  ultimate  numerical  evaluation 
should  be  kept  in  mind  when  the  mathematical  analysis  is  carried  out.  It 
is  a  definite  advantage  when  the  computations  required  by  a  theory  may 
be  readily  performed.  It  will  be  seen  that  when  hydrodynamical  problems 
are  attacked  by  the  operator  method,  the  formulas  which  occur  contain  a 
number  of  standard  functions  which  are  the  same  in  all  problems  of  the 
same  type.  Consequently  these  functions  may  be  computed  and  tabulated 
once  and  for  all.  When  this  has  been  done,  the  computations  required 
for  any  specific  problem  are  much  shorter  and  simpler. 


The  Method  of  Operators  in  the  Theory 
of  Compressible  Fluids 

The  idea  of  the  method  of  operators  is  now  explained  briefly. 

It  has  already  been  stated  that  the  partial  differential  equations 
of  compressibility  are  nonlinear  and  that  a  mathematical  theory  of  such 
equations  does  not  yet  exist.  Fortunately,  it  is  possible  to  linearize 
the  equations  of  motion  of  a  compressible  fluid  by  a  change  of  variables. 
Therefore,  hydrodynamical  problems  are  not  studied  in  the  physical  plane 
where  the  flow  occurs  but  must  be  investigated  in  some  auxiliary  plane 
where  the  differential  equations  are  easier  to  handle.  This  entails 
considerable  distortion  of  the  flow  pattern  and  a  consequent  loss  of 
intuitive  appeal. 

The  process  of  linearization  is  conveniently  performed  in  two  stages. 
In  the  first  place  it  is  noted  that,  in  general,  in  a  sufficiently  small 
portion  of  the  flow  region,  a  particle  may  be  just  as  unambiguously 
specified  by  its  velocity  vector  as  by  its  position,  although  this  will 

3m- 

The  main  mathematical  idea  of  the  operator  method  consists  in 
employing  certain  integral  operators  which  transform  functions  of  one 
complex  (or  real)  variable  into  "classes"  of  complex  solutions  of  given 
linear  partial  differential  equations  and  preserves  many  properties  of 
the  class  of  functions  to  which  the  operator  is  applied.  Thus,  in  the 
case  of  equations  of  elliptic  type,  the  application  of  the  operators 
results  in  the  derivation  of  complex  solutions  of  the  given  solution’s 
equation  which  in  many  respects  behave  like  analytic  functions  of  a 
complex  variable  and  help  in  the  investigation  of  real  solutions  in  a 
manner  which  bears  a  close  analogy  to  the  role  of  analytic  functions  of 
a  complex  variable  for  real  harmonic  functions  in  two  variables. 
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not  generally  be  true  over  the  flow  as  a  whole  .  The  first  change  of 
variables  is,  therefore,  from  the  physical  plane  with  coordinates  x,y 
to  the  so-called  velocity  or  holograph  plane  in  which  the  polar  coordi 
nates  q,0  are  used,  q  being  the  speed  of  the  particle,  and  ,  e 
angle  between  the  direction  of  motion  of  the  particle  and  the  positive 
x-axis.  This  transformation  obviously  introduces  distortion  (see  iol 
lowing  figs.)  and  it  is  by  no  means  bi— unique. 


y 


Motion  in  the  physical 
plane. 


The  image  in  the  hodograph 
plane  of  the  flow  indicated 
in  figure  at  left. 


since  it  usually  happens  that  the  velocity  vectors  are  equal  at  different 
points  of  the  flow  region.  A  further  distortion  is  now  produced  by 
plotting  log  q  end  0  as  Cartesian  coordinates  (logarithmic  plane). 
Finally,  the  substitution 


A  =  A(q)  =  £  arc  tan 


[h(M2  -  l)1//2]  -  arc  tan  [(m2  -  l)^  ] 


(see  fig.  l)  where 


k  -  iy/2 

k  +  1/ 


^Xn  order  to  make  the  correspondence  between  the  flow  pattern  in 
the  physical  plane  and  its  image  in  the  hodograph  plane  1-to-l,  it  is 
necessary  in  this  case  to  locate  the  hodograph  on  a  Riemann  surface  with 
the  proper  number  of  sheets. 
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has  the  effect  of  changing  the  scale  along  the  log  q-axis  (pseudo- 
logarithmic  plane).  This  introduces  an  additional  distortion  but 
succeeds  in  reducing  the  differential  equation  for  the  potential  func¬ 
tion  to  the  form 


where 


0aa  ~  0ee  -  -  0 


8  (M=-1)3/2 


(7) 


(8) 


The  stream  function  satisfies  the  differential  equation 

*aa  -  +  '•“l+A  '  0  (9) 

In  the  subsonic  case_,  a  similar  procedure  leads  to  the  differential 
equation 


*xx  +  tee  +  ^x  =  0 


do) 


(see  equation  (47)  of  reference  l),  where  X  =  X(q)  is  a  function  of 
the  speed  alone the  form  of  which  is  very  similar  to  that  of  A(q)  in 
equation  (5).  The  operator  method  of  finding  subsonic  compressible  flow 
patterns  is  now  easily  described.  It  is  based  on  making  suitable  use  of 
the  following  formula  for  generating  solutions  of  equation  (10): 


+(M)  =  ^  p(g) 


P(g)  =  H(2X) 


6(0 


y  qW(2x) 

^ ^  p2n  y.  t 

n=0  ^  11  • 


pc 


pc 


u  0 


>  (1 


nr-1 


S(?n)4tn 


where  £  =  X  —  i0.  (See  reference  1_,  formula^  (66).)  In  this  formula 
g  represents  an  arbitrary  analytic  function  of  a  complex  variable  £ 

^Note  that  in  formula  (66)  of  reference  1^  there  is  a  misprint:  the 

_  ,  (2n)!  ,  n  -  ,  (2n)! 

factor  - -  should  be  replaced  by 


pn 


nj 


2^n  n! 
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which  is  regular  at  t,  =  0;  whereas  H  and.  the  Q  n  ’s  are  certain 
functions  of  X  only  and  are  the  same  for  all  stream  functions.  They 
may  he  computed  once  and  for  all  and  the  results  used  in  all  future 
hydrodynamical  problems.  (See  reference  4.) 

The  foregoing  formula  is  employed  in  the  following  manner.  If  any 
analytic  function  g  of  a  complex  variable  (regular  at  £  =  0)  is 
substituted  into  formula  (ll),  the  result  is  a  stream  function  of  a 
subsonic  flow  pattern.  On  the  other  hand,  it  is  well  known  that  an 
analytic  function  g(log  q  -  ±6)  may  be  regarded  as  defining  a  flow 
pattern  of  an  incompressible  fluid.  Consequently  formula  (ll)  may  be 
interpreted  as  an  operator  which  distorts  an  incompressible-flow  pattern 
into  a  compressible  one.  (it  may  be  said  that  the  operator  acts  as  a 
distorting  mirror  which  reflects  an  incompressible-flow  pattern  as  a 
subsonic  compressible-flow  pattern.)  The  operator  formula  is  of  such  a 
form  that  a  knowledge  of  the  incompressible— flow  pattern  used  leads 
directly  to  both  qualitative  and  quantitative  knowledge  of  the  behavior 
of  the  stream  function  of  the  corresponding  compressible  flow. 

An  obvious  objection  is  that,  whereas  any  incompressible  flow  will 
lead  to  some  compressible  flow  when  a  specific  obstacle  is  prescribed, 
there  seems  to  be  no  indication  of  how  to  find  precisely  the  g  to  be 
substituted  into  the  operator.  A  procedure  to  solve  this  problem 
approximately  can  be  developed  (see  reference  7)*  hut  this  is  not  con¬ 
sidered  herein.  However ,  it  is  easy  to  see  why  the  same  analytic 
function  g(£)  which  solves  the  corresponding  incompressible  problem 
can  be  expected  to  yield  a  first  approximation  to  the  compressible 
problem  in  the  case  of  low  speeds  (i.e.,  small  Mach  number);  for  in  such 
a  case  the  solution  (formula  (ll))  is  easily  shown  to  yield  a  flow,  very 
similar  to  the  incompressible  flow,  which  agrees  with  the  physical  state 
of  affairs. 

The  operator  method  as  just  described  involves  the  following 
elements : 

(1)  It  is  necessary  to  have  considerable  knowledge  about 
incompressible— flow  patterns ,  in  particular  about  their  images  in  the 
holograph  and  logarithmic  planes . 

(2)  The  operator  formula  has  to  be  derived,  the  functions  H 

and  have  to  be  determined,  and  rules  must  be  developed  for 

interpreting  properties  of  incompressible— fluid  motions  as  properties 
of  the  compressible  flows  given  by  the  operator. 

(3)  After  a  solution  of  a  problem  on  compressibility  has  been  found 
in  terms  of  the  variables  (X,0)  (i.e.,  in  the  so-called  pseudo- 
logarithmic  plane),  means  must  be  furnished  for  determining  the  actual 
flow  in  the  physical  plane. 

It  is  obvious  that  one  reason  why  this  method  can  be  so  powerful 
fn  the  study  of  subsonic  compressible  flows  is  that  the  theory  of 
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incompressible  flows  (i.e.,  the  theory  of  analytic  functions  of  a  complex 
variable)  is  so  well  known. 

The  supersonic  case  is  attacked  by  perfectly  analogous  means.  As 
has  been  indicated  by  the  author,  there  exist  various  operators  for 
transforming  solutions  of  one  equation  into  solutions  of  another  equa¬ 
tion;  and  often  for  different  purposes  it  is  convenient  to  use  different 
operators  or  different  forms  of  the  same  operator.  In  the  present 
report  two  operators  which  transform  solutions  of  the  relatively  simple 
hyperbolic  equation 


=  0  (12) 

Sa2  Be2 


(i.e.,  a  pair  of  differentiable  functions 


g(?)  and  3a(  r|) 

where  £  =  (A  +  0)  and  r\  =  (A  —0) 


(13) 


of  one  variable)  into  solutions  of  the  compressibility  equations  are 
derived. 

As  soon  as  the  function  \|r(A,@)  satisfying  equation  (9)  is  obtained, 
the  corresponding  flow  in  the  physical  plane  can  be  obtained  exactly  in 
the  same  manner  as  in  the  subsonic  case.  (See  reference  1,  section  14.) 

In  contrast  to  the  subsonic  case,  flows  which  in  the  logarithmic 
plane  satisfy  equation  (12),  A  being  some  function  of  the  speed,  have 

not,  to  the  knowledge  of  the  author  of  this  report,  been  studied. 
Relatively  little  is  known  concerning  how  to  choose  g  and  h  in  order 
to  obtain  in  the  physical  plane  a  flow  around  a  prescribed  obstacle, 
although  there  is  reason  to  expect  that  this  investigation  should  be 
simpler  than  in  the  subsonic  case.  However,  this  question  is  not  con¬ 
sidered  herein,  but  is  postponed  to  a  future  publication. 

The  following  sections  of  this  report  are  devoted  exclusively  to 
the  determination  of  two  different  operators  which  transform  differ¬ 
entiable  functions  of  one  real  variable  into  solutions  of  equation  (9). 

^In  reference  8  Coburn  has  indicated  the  extension  of  the  Chaplygin— 
Karman— Tsien  method  to  the  case  of  supersonic  flows.  The  stream  functions 
of  such  flows  satisfy  equation  (9).  Coburn  indicated  one  characteristic 
property  of  these  flows.  On  the  other  hand,  the  general  properties  and 
the  form  of  flows  (in  the  physical  plane)  which  are  of  interest  for  the 
purpose  of  the  present  approach  have  not  been  investigated  in  Coburn1 s 
paper . 
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Operator  Obtained  by  the  Use  of  Riemann's  Function 

One  method  of  generating  stream  functions  of  supersonic  flows  from 
expression  (13)  is  by  the  use  of  Riemann's  function. 

As  indicated  in  the  preceding  section,  the  stream  function  ty(A ,0) 
satisfies  equation  (9) .  If,  instead  of  A  and  0,  the  variables  | 
and  n  defined  in  expression  (13)  are  introduced,  and,  instead  of  \|r, 
the  "reduced  stream  function" 


i*  =  tyexp 


(14) 


(where  a  is  an  arbitrary  constant)  is  considered,  then  a  formal  compu¬ 
tation  shows  that  i|r*  satisfies  the  equation 

ir*  +  F-J?  +  =  0  (15) 

where 


m  2  _  =  k  +  1  f~5(k  +  l)  12k  6k  -  lk 

1  1  d(2A)  6k  [  b6  ^  b2 


+  (kk  +  8)  -  (3k  -  1)B2J  (16) 

where 

B2  =  M2  -  1  (17) 

Remark.—  Obviously  it  is  allowed  to  take  for  the  lower  limit  a 
of  the  integral  in  equation  (lk)  an  arbitrary  point,  since  replacing  a 
by  a*,  a*  ^  a,  means  only  that  is  multiplied  by  a  constant,  . 

Nb(t)  dT .  The  most  natural  choice  is  to  take  for  a 

the  same  value  as  in  the  subsonic  case.  (See  formula  (ill)  of  reference  1.) 

As  is  explained  in  more  detail  in  appendix  C,  it  is  possible  to 
extend  the  variable  A  to  the  complex  values  1  -  X  +  iA.  Then  in  the 
subsonic  case  the  variable  l  assumes  values  on  the  negative  real 
axis  X  and  on  the  positive  imaginary  axis  A  in  the  supersonic  case. 

The  integration  can  be  then  taken  from  t  =  -00.  Since  N(0)  =  °°,  the 
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integration  must  "be  performed  along  a  path  which  avoids  7  =  0,  for 
example,  along  the  real  axis  from  —00  to  X <  0;  then  along  a 

circle  with  radius  |  Xy  |  and  center  at  2  =  0;  and  finally  along  a 
segment  of  the  positive  imaginary  axis  from  r  =  i|XjJ  to  t  =  jLA. 

Let  ,  p;  ?0,t]0)  be  the  Niemann’s  function  of  equation  (15) . 

(See,  e.g.,  pp.  311—317  of  reference  9-)  Then,  according  to  a  classical 
result  (p.  316,  formula  (7’)  of  reference  9), 

**(M)  =t*(i0^o)  E(Mi60,T|o) 


+  **y(^Y)  dY 


^*x(Y->rlo)  E(X-,r|o^o> %)  dY 


where  i|r*j  =  r-L,  and  so  forth. 
oX 


Since  for  i|r*(l0,  t)q)  can  be  chosen  an  arbitrarily  prescribed 
constant  A  and  for  ->|r*(£0,Y)  and  t)0)  can  be  chosen  two  arbi¬ 

trary  continuous  functions  h(X)  and  g(Y)  ^provided  that  h(iQ)  =  S(T10)^ 
the  formula 


+*(6,1)  =  AN(I,ti;|0,ti0)  +  h(Y)E(|0,Yj|0,T,0)  dY 


g(X)E(X,r)o;|o,q0)  dX 


represents  an  integral  operator  which  transforms  two  arbitrary  continuous 
functions  of  a  real  variable  into  a  solution  of  equation  (15). 
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(M) 


(S0»  'np) _ _ ^,T]0^ _ x 

The  rectangle  within  which  the  values  of  **(M)  are  determined  by 
prescribing  its  values  along  the  line  segments:  Y  =  ~  ^  — 

and  X  =  S0,  Tl0  —  Y  —  T>- 


As  is  well  known,  the  Riemann  function  of  equation  (15)  can  he 
written  in  the  form 


00 

R(£5TU!0,T10)  =  ]>~  (-l)n  R<n)( 6,ijjl0,il0) 


where 


:(o)(6,tli60,Tl  )  =  1 


H  p 

R(i)(|,t1;|0,t1o)  =  f±(x  +  y)  ax  aY 

Li^o  <1% 


I  flT| 


(2)(I^;I0,ti0)  =  FX(X  +  dX  dY 


i  m 


F,(X  *  T)E(“-X)(I,T;S0,ti0)  4X  dY 


0% 


(See  reference  10.) 
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Except  for  the  fact  that  Fj_(0)  =  00,  which  causes  certain  complica— 
4  tions,  the  evaluation  of  the  foregoing  integrals  does  not  represent  any 

essential  difficulty.  The  tabulation  of  E  =  R(£,t];£0,t)0)  is  rather 

complicated  since  R  is  a  function  of  four  variables. 

In  appendix  A  the  second  term  R^  of  the  expansion  (equation  (20)) 
is  computed  and  tabulated. 

Remark.—  In  reference  1,  formula  (119)^,  the  expression  for  a  funda¬ 
mental  solution  of  equation  (121 )  has  been  derived.  The  coefficient  a 
of  the  logarithmic  term  is  given  by  equation  (123).  It  has  exactly  the 
same  structure  as  the  expression  for  R(|, rj;  |Q, tj0),  except  for  the  fact 

that  the  arguments  of  a  are  complex  quantities,  whereas  £jT)J£0<,‘r)0 

are  real.  In  reference  11,  Kraft  has  elaborated  a  working  procedure  for 
the  evaluation  of  a  which  can  be  used  for  actual  computation  of  R. 

It  becomes  in  this  case  simpler  since  the  arguments  are  real. 


Integral  Operator  of  Type  Given  by  Formula  (ll) 

In  analogy  to  the  subsonic  case,  where  formula  (ll)  yields  a  repre— , 
sentation  for  solutions  of  equation  (10)  in  terms  of  an  arbitrary  analytic 
function  of  one  complex  variable,  there  is  derived  in  this  section  a 
representation  distinct  from  the  one  explained  in  the  preceding  section 
for  solutions  of  equation  (9)  (the  equation  for  the  stream  function  in 
the  supersonic  case)  in  terms  of  two  arbitrary  differentiable  functions 
of  a  real  variable. 

Obviously,  instead  of  equation  (9),  equation  (15)  for  the  reduced 
stream  function  -\|f*  may  be  considered.  The  desired  representation  can 
be  obtained  as  a  consequence  of  the  following  theorem. 

Theorem^.—  Let  E^n)(T)  denote  functions  defined  by  the  following 
recurrence  relations: 


E 


(i)(T)  .  _ 


Ot 


fi(ti)  ari 


V  (22) 


I  f  _ 

e(»+i>(t)  .  -  [E<n>  (Ti)  +  Fi(Ti)E<n)(Ti) 


dT 


^Note  that  in  equation  (119)  and  line  10  of  page  46  of  reference  1, 
by  a  misprint,  a  capital  A  instead  of  a,  lower  case,  is  used. 

®This  theorem  has  been  indicated  in  reference  3;  see  chapters  4 
and  5  of  this  reference.  See  also  reference  12. 
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Here  F-^t)  is  a  given  function,  the  derivatives  of  which  satisfy  the 
inequalities 


dnF!(T ) 


dT 


n 


< 


(n  +  1)!M 


(To-  e  -T) 


n+2- 


n  =  0,1,2, 3, 


(23) 


where  e(>0),  M  =  M(e),  and  tq  are  conveniently  chosen  constants. 
Then  the  series 

00 

=  f(|)  +  Y_  E(n)^  +  ^)f(n)(0  (24) 

n=l 

and 

00 

=  e( ti)  +  ^_  E^(£  +  t)g^ntri)  (25) 

n=l 


where  f  and  g  are  two  arbitrary  differentiable  functions  of  a  real 
variable  and 


fk+o-ta) 


fM(li)  ilj 


> 


(26) 


g  Cn+1]  (tj) 


IT] 

Jo 


g^(T1l)  dT]l 


» 


I 


will  converge  in  the  intersection  of  the  domains 


-1  < 


To  “  1  "  1 


<  1 


and 


-1  < 


<  1 


(27) 


I 

and  will  represent  there  solutions  of  the  equation 


H(\|r*)  =  +  F1(|  +  q)t*  =  0 


(28) 
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The  proof  of  this  theorem  is  given  in  appendix  B.  As  has  been  shown  in 
t  reference  13  (see  equation  (45) ),  F^  is  given  by  equation  (l 6) 9 

where  B^  and  2 A  are  connected  by  relations  (5)  and  (17)*  Thus 

the  question  of  the  representation  of  solutions  of  equation  (15)  is 
reduced  to  the  investigation  of  the  behavior  of  the  function  F1(t)j  in 

particular  to  determining  whether  its  derivatives  satisfy  the  inequali¬ 
ties  (23).  This  investigation  has  been  carried  out  in  appendix  C.  If 
in  the  Ay6— plane  the  origin  is  shifted ^  thus  introducing  instead  of  2A 
the  new  variable 


and  therefore  in  particular  in  the  circle  of  radius  tq  —  e  = 

in  (I  --  1^  —  ^  e  >  0.  (See  fig.  2.)  If  M  =  M(c)  denotes  the 
maximum  of  F^(2A)  in  this  circle^  then  by  classical  results  of  the 
theory  of  functions  of  a  complex  variable^  F^  satisfies  the  inequali¬ 
ties  (23),  where  2Af  is  substituted  for  t#  Since  e  can  be  chosen 
arbitrarily  small ^  the  series  (equations  (24)  and  (25))  formed  for  equa¬ 
tion  (31)  converge  in  the  domain  indicated  in  figure  3* 

On  the  other  hand,  in  many  instances  it  is  necessary  to  have  repre¬ 
sentations  (equations  (24)  and  (25))  of  the  solutions  which  are  valid  in 
the  whole  strip^ 


—00  <  6  < 


00 


0  <  2A  <  2n(±  -  l) 


(33) 
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One  possibility  of  obtaining  such  formulas  consists  in  approxi¬ 
mating  F2(2A* )  by  a  sequence  of  functions  F2(m)(2A* each  of  which 

is  an  entire  function  of  A  and  therefore  (by  a  classical  result  of  the 
theory  of  analytic  functions  of  a  complex  variable)  satisfies  the 
inequalities  (23)  and  which  is  chosen  in  such  a  manner  that 


lim  F2W(2A')  =  F2(2At) 


in  the  interval 


0  <  2A  <  2(j-  -  ljjt 

(Concerning  another  method  of  representing  Y^z ' ,  k  =  1,2,  see  (27),  (25) 

in  the  domain  (33),  and  reference  2,  section  5.)  Let  (|  * , rj* )  denote 

the  solutions  of  the  equation 


where 


+  F2(m)(S'  +  T|’ )**<“>  -  0,  m  -  1,2, 


-s(0)  <  S’  <  s(0! 


such  that  \]f*(in)  (|  «^0)  =  ty*(| 1 ,0)  and  \(r*(m)  (0, -rj * ) 
_n(0)  <  tj *  <  ^(0).  As  is  shovn  in  appendix  E, 


=  t*(0jT1‘), 


lim  (| '  ,r)‘ )  =  ^*(| f  j T|* ) 


in  the  interval 


(0)  ,  (0)  (0)  ,  (0) 

-5  <  |  <  5  ,  — n  <  -n*  <  n 


A  sequence  of  functions  F, 
manner .  Let 


^m)(2A‘)  can  be  obtained  in  the  following 


F2(2A* )  =  aQ  +  cc^A*  +  <x2A'  +  .  .  . 


(39) 
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be  the  series  development  of  Eg  around  the  point  2A»  =  0 
(l.e.,  2A=(i-l)ny  The  series  J  o^A*11  obviously  converges 


in 


n=0 


the  circle  indicated  in  figure  2.  The  series 


2W(2A')  = 


n=0 


r(l  +  n/m) 


-  A* 


n 


(40) 


J°r  finite  positive  number  m,  an  entire  function  of  A' 

and  therefore  the  series  (equations  (24)  and  (25))  obtained  for  equa¬ 
tion  (36)  with  F2(m)  given  by  equation  (40)  will  converge,  for 

any  m  <  »  in  the  strip  (33).  On  the  other  hand,  according  to  a 
classical  theorem  of  the  theory  of  functions  (see  reference  14),  the 
relation  (equation  (34))  holds  in  this  case,  so  that  for  the  functions 

tnifl  obtained  in  the  foregoing  manner  the  relation  (equation  (37)) 

holds  in  the  vhole  interval  (38). 

Naturally  it  is  also  possible  to  approximate  F2(2At)  by  using 

different  procedures.  Some  of  them  are  convenient  for  some  special 
purposes,  in  particular  if  a  problem  requires  the  representation  of 
solutions  in  certain  subdomains  of  the  strip  (33). 


SUMMARY  REMARKS 


The  present  paper  employs  the  holograph  method  for  generating  flow 
patterns  of  supersonic  flows.  The  essential  feature  of  this  method  is 
that  the  equations  are  linear  and  therefore  the  principle  of  superposition 
of  solutions  holds.  If,  therefore,  *v(A,e),  V  =  1,2,3,  .  .  f  Jepre^ 

eents  a  set  of  particular  solutions  of  equation  (7),  and  Av,  Arbitrary 
constants^  a ny  linear  combination 


N 


VvK0) 

V=1 


is  also  a  solution  of  equation  (7).  By  varying  the  constants  Av, 

IT  1,t>  A  \l’ *  flOWS  around  different  shapes  can  be  obtained.  On  the 
other  hand  it  is  often  necessary  to  determine  constants  Av  to  yield 
a  flow  which  approximates  that  about  a  prescribed  boundary  curve  the 
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equation  of  which  is,  say,  F(*,y)  =  0.  In  this  case,  the  constate  A 
can  "be  determined  in  such  a  manner  that 


M(A)  e 


1 

l; 


F[x(A,e;Av),y(A,e;Av)]  |  2  ae 


■under  t-he  condi'fcion 


=  0 

Here  x  =  x(A,e;Av)  and  y  =  y(A,@;Av)  represent  functions  of  A  and  6 
defined  in  the  INTEODUCTION  and  corresponding  to  the  stream  function 
if  =  y  Avifv(A,6)  in  the  A, 0-plane. 

«.  >•  -n  T-  rf  "  *“  s,s:s: 

=£vS:=: i 

zero  vhen  the  number  W  of  functions 

sidered  as  a  condition  for  the  existence  Ox  a  solution. 

On  tie  other  hand,  in  many  instances,  only  solutions  of  the  problem 

£  “ordered  hut  also  the  condition  that  theae  shoe*  conditions  are 
satisfied  along  unknown  characteristics. 


Brown  University 

Providence,  B.I., 


November  l8,  19^6 
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APPENDIX  A 

t 

EVALUATION  AND  TABULATION  OF  THE  FUNCTION 


I 


I 


The  second  term  in  the  series  for  the  Riemann 

function  is  evaluated  in  this  section. 

F1  is  a  function  of  the  variahle  A  =  ^  +  11  and  therefore 


Fx(x  +  T) )  dx  =  r(1)(!  +  ti)  - r(l)(i0  +  h) 


fir) 

dY 

%  U 


05 


OH 


F(X  +  Y)  dX 


[r(l)(i  +  Y)  -  r(1)(iD  +  Y) 


dY 


(K 


=  r(2)U  +  T|)  -  r(2)(?  +  1J0)  -  r(2)(io  +  n) 

+  rt2)(^o  +  \)  ^ 

where 

1ti 

f1(t)  dr  =  r^(tx)  -  r(^1^(to)  (42)' 

JtD 

an& 

r(1)(T)  dr  =  r^2^)  -  r(2)(t0)  (43) 

^o 
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The  fact  that  F(0)  is  infinite  causes  a  certain  amount  of  inconvenience 

in  the  tabulation  of  the  functions  f  and  r(2t  According  to  equa¬ 
tion  (l6),  for  k  =  1.4, 

F  (b)  =  hhl  +  _  0.51  -  0.12B2  (44) 

1  B6  B^  B2 

where 

A  =  2  \[6  arc  tan  —  2  arc  tan  B  (45) 

\|6 

and  B2  is  given  by  equation  (17),  M  being  the  Mach  number. 

In  table  I  the  values  of  M,  B,  and  F1  have  been  tabulated?  in 
figure  4,  F^_  is  given  as  a  function  of  M,  B,  and  2A.  For 
2A  =  O.3686  the  function  F^_  vanishes,  that  is,  F]_  [b(0.3686)]  =  0; 

this  point  has  been  chosen  as  the  initial  point  for  the  integration, 
since,  as  has  been  pointed  out  earlier,  for  the  most  natural  choice 
2A  =  0,  becomes  infinite. 

The  values  of  r(1)(2A)  and  f^2^(2A)  for  (O.3686)  = 

r(2) (0.3686)  =  0  are  given  in  table  II.  The  interval  of  tabulation 

of  2A  is  one  four-hundredth.  In  order  to  obtain  T^2^(2A)  for  inter¬ 
mediate  values,  the  interpolation  formula 

r(2)(2A)  =  r^2^(2AQ)  +  pD-j^  +  p2D2  (46) 

can  be  used,  where 

r(2)(2A)  =  r(2)(2AD  +  ph)  =  r(2)(2A0)  +  ph  +  |p  (p  -  l)h2 

=  r(2)(2AQ)  +  pD^  +  p2D2  (47) 

The  quantity  h  is  the  interval  of  tabulation  at  2A0<> 
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APPENDIX  B 

PROOF  OF  THEOREM  GIVEN  BY  EQUATIONS  (22)  TO  (28) 


The  proof  is  carried  out  at  first  under  the  assumption  that  the 
series^  (equations  (24)  and  (25))  as  well  as  those  for  their  first 
derivatives  converge  uniformly.  (The  validity  of  this  assumption  is 
proved  in  this  appendix.)  A  formal  computation  yields 


V^1)  ,  t  =  fE^1)  +  )  f M 

6  V  V 

n=l 


It],  +  E^n+1) 


(48) 


00 

F^1)  =  fF1  +  ^  fCn]p1E(n, 


n=l 


that  is,, 


(49) 


OO 

,  +  F-^1)  =  f  E^1)  t  +  fJ  +  ftn]J~E^n 

L  J  n=l  L 

+  E^n+1^,  +  F-^E^J 


tW 


(50) 


where  is  a  function  of  one  variable  t  =  2A1 . 

Assuming  that  the  functions  are  also  functions  of  one 

variable  ry 


E^n)  t  =  eW 


E^n^ . »  •  =  E^n) 


> 


(51) 


s  H  "  tt 

which  yields  the  foregoing  relations 9  if  the  assumption  is  made  that 

E^n)(a)  =  0 


In  this  appendix  are  replaced  by 


(52) 
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It  remains  therefore  to  prove  that  the  series  and  both  its  formal 
(i.e.j  term-by— term)  derivatives  converge  uniformly. 

Notation.—  If  in  the  interval  I  there  hold  for  the  functions  A(t  ) 
and  2(t)  and  all  their  derivatives  the  inequalities: 


|a(t)|  <  $(t)  and 


dnA(T )  <  d°%(T) 

dTn  dTn 


(53) 


A(t)  will  be  denoted  as  a  dominant  of  A(t),  which  fact  is  symbolized 
by  writing  A «  A  or 


A  »  A,  re  I 


(5^) 


If  e(1Ht)  is  given  by 


e(D  (t) 


At 


\(Tl)  dTl 
do 


(55) 


where 


Fx(t)  »  F1(t ) 


(56) 


then 


^(^(t) 


At 


a 


Fl(Tl)  dTl 


0T 


!1(t1)  dr  <  E^fr),  Tel  (57) 


t/a 


and  also 


Thus 


dnE(1)(r) 

dTn 


^  dnE(1) (t) 

~  „  9 

dTn 


T€I 


(58) 


e(1)(t)  «  e(1)(t),  Tel 


(59) 
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Suppose  now  that 


nr 


E(n+1)(r)  = 


lf(n)  ( T 

T 


,T1(Tl)*il(Tl)i,,l|’l)  4t1  +  S<“>(t)  (60) 


where 


E 


^n)(T)  «  E^t),  0  «  §(n)  (t)  }  Tel 


(61) 


Then  it  follows  Immediately  that 


E(n+l)(T)j  <  E(n+l)(T) 


dE(n+l)(T) 


dT 


< 


dE(n+l)  (t) 
It 


(62) 


and^  hy  considering  the  corresponding  derivatives  of  [S^n^TT  + 
in  comparison  with  [E^n)TT  +  F^11)],  it  follows  that 


E(n+l)(T)  «E(n+l)(T)^  Tej 


(63) 


which  completes  the  proof  hy  induction. 
Now  hy  expression  (23), 

fi(t)  «  F^t)  = 


M 


(To~e-T) 


2n 


,  Tfl 


(64) 


If  H^t)  is  given  hy 


jj(n)  =  n  o 


(To  "  e  “  a) 


n 


(65) 
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where  the  C  's  are  some  conveniently  chosen  positive  constants  (to  he 
n  /v» .  . 

determined  later) ,  then  an  explicit  expression  may  he  obtained  for  E(t), 
namely., 

:(■](,)  = — —  (66) 
(To  -  6  -  T)“ 


In  order  to  express  Cn+^  in  terms  of  Cn>  the  right-hand  side  of 
equation  (66)  is  substituted  into  equation  (60)  to  give 


g(n+l)  (x)  =cnMTQ: 


IT 

n(n+  1)  +  CnMr02 

/  ,  .n+2 

\j  a 

L(to  -  e  -  T)  J 

dT  + 


Cn+1MV 


Cn+lMTo 


where 


(67) 


Cn+1  =  Cn  ln  + 


Mt. 


n  +  1 


(68) 


Thus  for  sufficiently  large  values  of  n 

cn+l  =  Cn(n  +  D  (69) 

and  therefore  for  every  value  of  M 

Cn+1  ^  M*(n  +  1)!  (70) 

where  M*  is  a  conveniently  chosen  constant.  The  function  f (?)  is 
assumed  to  he  differentiable  and  therefore  there  exists  a  constant, 
say,  Mq,  such  that 


» 


1 


4 


|  f  (i* )  |  <  M2  for  1^*1  =  h 


Then 


fM(l’)| 


<^,n 

=  n! 


(71) 


* 


for  | ? ’ |  —  ^2 


(72) 
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APPENDIX  C 

INVESTIGATION  OP  THE  FUNCTION  F1(2A) 


This  section  is  devoted  to  the  investigation  of  the  behavior  of 
the  function  F1(2A)  in  order  to  determine  where  its  derivatives  satisfy 

the  inequalities  (64). 


1  /2 

If  M  varies  from  1  to  B  =  (m2  —  l)  '  varies  from  0  t 


A  =  i  arc  tan  (hB)  —  arc  tan  B 
h 


from  0  to  —  1^.  Therefore  it  is  necessary  to  investigate  the 
behavior  of  F^PA)  on  the  interval 

s  - 1)«]  w 

Since  F1  (given  by  equation  (l6))  is  a  rational  function  of  B,  it  is 
sufficient  to  investigate  the  function  B(2X). 

In  order  to  carry  out  this  investigation  it  is  convenient  to  continue 
the  functions  to  complex  values,  so  that  it  will  be  possible  to  use 
methods  of  the  theory  of  analytic  functions  of  a  complex  variable. 

Instead  of  equation  (5),  the  complex  function 


2(3)  =  X  +  iA  =  \ 


H-iiCg 


(78) 


of  the  complex  variable  P  =  T  +  IB  will  he  considered  ^p  =  T,  0  <  T  <  1 

When  p  is  real  it  is  seen  that,  if  the  principal  values  of  the 
logarithmic  terms  are  chosen,  then  l  is  also  real  and  the  foregoing 
equation  reduces  to: 
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1  -  T  1  .  h-1  +  T 

log  -  +  —  log  - 

>  1  +  T  h  h-1  -  T 


,  A  =  0 


in  agreement  with  the  earlier  definition  of  X.  On  the  other  hand,  if 

3  =  iB,  B  >  0 

B  is  purely  imaginary;  then,  if  the  principal  values  of  the  logarithmic 
terms  are  taken  once  again,  there  results: 


arc  tan  B  +  —  arc  tan 
h 


(hB)] 


which  is  also  in  agreement  with  the  previous  definition.  The  question 
which  has  to  he  investigated  is  to  determine  the  domain  of  regularity 
of  3(2),  inverse  to  equation  (77). 

A  classical  theorem  concerning  the  inversion  of  analytic  functions 
states  that,  if  w(z)  is  a  function  of  the  complex  variable  z  which 
is  regular  at  z  =  zQ  and  has  a  nonzero  derivative  there,  then,  in  some 

sufficiently  small  neighborhood  of  the  point  w  =  wQ  =  w(zQ),  it  is 

possible  to  invert  the  function  w(z),  that  is,  to  express  z  as  a 
regular  function  of  w.  (See  reference  15,  p.jA2.)  At  every  point 
except  3  =  ±1  and  3  =  ih,  the  function  2(3)  is  regular,  so  that 
the  zero  points  of  the  derivative  d2/d3  have  to  be  determined. 

A  formal  computation  yields 


1  _  1  1  f _ 

1-3  1  +  *3  +  h  l  h~ 1 


h~ x  +  3 


(l  -  h2)^2 
(l  -  3^)(l  -  h2^2 


Obviously, 


—  -  0  for  3=0  and  3 
d3 


The  discussion  would  become  quite  complicated  if  alLthe  possible 
branches  of  the  function  2(3)  and  of  the  inverse  3(2)  were  to  be 
considered,  because  each  summand  of  equation  (8l)  is  an  infinitely 
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many-valued  function  of  P,  having  logarithmic  singularities  at  +1 

and  ±h-1,  respectively.  It  is  therefore  convenient  to  render  2(P) 
single— valued  hy  the  artifice  of  starting  at  P  =  0  with  the  principal 

values  of  log(l  -  p),  log(l  +  P),  logjr1  +  p)  and.  logClT1  -  p) 

(i.e.,  with  the  real  values  of  the  logarithms)  and  slitting  the  complex 
(1-plane  from  1  to  +°°  and  from  —1  to  —  °°  along  the  real  axis. 

Since  h-1  >  1,  it  is  clear  that  the  slit  plane  contains  no  branch 
points  and  therefore  2  is  uniquely  defined  hy  beginning  with  the 
aforementioned  determination  of  2( 0) . 

A# 

Now,  a  direct  investigation  shows  that  the  |3-plane,  slit  in  the 
manner  just  described,  is  mapped  into  the  band 

—oo  <  A.  <  <»,  — (h— ■*"  —  l)«  <  A  <  (h  ^  —  l)jt  (83) 

However,  the  mapping  is  not  hi — unique,  for  at  the  branch  point  p  —  0 
the  mapping  is  easily  seen  to  be  3  to  1.  (See  appendix  D.)  Thus,  the 
simply  covered  slit  'p-plane  is  mapped^into  the  simply  covered  band  just 
defined.  The  only  singularities  of  P(2)  in  the  strip 

— 00  <  X  <  00,  0  <  A<  ( h— —  l)it  (84) 

are  the ^points  2=0  where  P(2)  has  a  branch  point  and  2  =  (h-1  -  l)*i 
where  p(Z)  has  a  pole. 
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APPENDIX  D 


THE  DEVELOPMENT  OF  p(Z)  AT  THE  ORIGIN  Z  =  0 


The  simplest  method  of  proving  that  1  =  0  is  a  branch  point  of  $(Z) 
of  third  order  is  to  determine  the  series  development  of  this  function 
at  Z  =  0.  In  order  to  obtain  this  development,  the  principal  branch 
of  Z(p),  namely. 


\  log  +  h  108  =  ijj-  arc  tan  (p/i)  +  ^  arc  tan  (Ph/i)j  (85) 


has  to  be  developed  in  a  power  series  to  give 

Z  =  -  ^  (1  _  h2)  -  ^  (l  -  hM  -  £  (1  _  h6) 
3  5  7 


1  -  h‘ 


l  -  *r3  ■  3  (l-h4)  p  3  (l  -  h6) 

'  5(l_h2r  +  Ml-h2) 


Introducing  a  new  variable  s, 


83  - - 3 -  2,  that  ie_,  s  =  ( - - - 

1  -  h2  V  1  -  h2  .  / 


yields  the  development 


A# 

s  as  CXP  +  C^P3  +  C^P^  +  .  . 


Since  Cj_,  as  is  easily  seen_,  is  not  zero,  it  follows  from  the 
general  theory  of  analytic  functions  that  this  series  can  he  inverted: 


p  =  d^s  +  d-28^  +  +  .  .  • 


(90) 
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In  particular,  for  h.  =  l/\|6, 

o±  =  1,  c3  =  0.2333,  ...  d-L  =  1,  d3  =  -O.2333  •  ..  .10  (91) 


£ 


10The  values  of  cn  and  dn,  n  =  1,  3,  .  .  . , 
puted  "by  the  author  and  will  he  published  elsewhere. 


25  have  been  com— 
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APPENDIX  E 


APPROXIMATIONS  OF  SOLUTIONS 

In  the  present  report,  solutions  of  equation  (15)  have  been  approxi¬ 
mated  at  several  places  "by  solutions  of  an  equation  in  which  the  coeffi¬ 
cient  F-j_  is  replaced  by  another  coefficient  which  approximates  it.  The 

legitimacy  of  such  a  procedure  is  established  by  the  following  theorem. 
Theorem.—  Let  i|/k)  ( l,T)) ,  t  =  1,2,  be  solutions  of  the  equations 

tft  M  +  f<kVk)  =  0  (92) 

where  the  f(k)’s  are  given  functions  of  |  and  T)  which  satisfy 
throughout  the  rectangle  0  <  |  <  a,  0  <  t]  <  b  the  inequalities 

|  f^  —  f^  |  <  e  and  |  f^k^  |  <  m  (93) 

€  and  m  being  positive  constants.  The  functions  '  are  also  to 
satisfy  the  conditions: 

+^(s,0)  =  ^2^(fi,0)  =  X1(|)  for  0<?<a  (94) 

and 

t(l)(0,T,)  =  t(2)(0,Ti)  =  XgU)  for  0<T)<b  (95) 

where  X^  and  X2  are  continuous  functions  prescribed  on  the  closed 

intervals  (0,a)  and  (0,b),  respectively,  subject  only  to  the  condi¬ 
tion  that  Xx(0)  =  X2(0).  Let  a  be  a  positive  number  such  that 
|  Xj  |  <  a,  |  X2  j  la,  Then  there  exists  a  positive  number  M,  depending 

only  on  e,  a,  m,  and  ab  (the  area  of  the  given  rectangle)  such  that 
for  fixed  values  of  a,  m,  and  ab,  M  approaches  a  finite  (positive) 
limit  as  e  approaches  zero,  and  such  that,  throughout  the  given  rectangle, 
the  following  inequality  holds: 


t(l)(l,T,)  -  ^2^(I,T))  I  <  €M  (96) 


Proof.—  According  to  the  general  theory  of  linear  hyperbolic  partial  differential  equations,  it 

is  known  that  the  functions  satisfying  expressions  (92),  (93),  and  (95)  may  he  expressed  in 

the  form  of  the  following  infinite  series,  which  is  valid  throughout  the  given  rectangle: 
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so  that,  according  to  expression  (93),  |  5(4,^)  |  <  e.  Now  set  k  equal  to  1  and  2,  respectively, 

in  equation  (97)  and  subtract  the  two  equations  so  obtained  (pairing  off  corresponding  terms  and 

ing  f*(-0  by  p(^)  +  S  The  result  obtained  is: 


+  T-lip  "t|p  g 
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nrirj  it  is  clear  how  a  similar  reduction  can  "be  effected  in  all  the 

succeeding  terms ,  so  that  no  sums  appear  in  any  of  the  integrands.  #  i 

Nov,  if  each  term  in  the  series  on  the  right-hand  side  is  replaced 
by  its  absolute  value,  an  upper  bound  is  obtained  for  the  absolute  value 
of  the  left-hand  side  of  equation  (100);  and  this  inequality  continues  i 

to  hold  if  the  functions  p,  8,  and  f  ^  are  replaced  by  their  respec-  j 

tive  upper  bounds  3^^  and  m.  (The  upper  hound  for  p  follows  [ 

directly  from  equation  (98)  and  the  definition  of  a.)  Therefore^  from 
equation  (100)  there  is  obtained  the  inequality:  !i 

(3!)2 


It  is  not  difficult  to  see  that  the  nth  term  of  the  right-hand  side 
of  equation  (102)  is  simply: 

3a  [(e  +  m)n  -  mj|nr)n 
(n!)2 


(103) 

f 


♦ 


(1)  _  t(2) 


<  3<X€|.T1  + 


3a(e2  +  2em)|2T]g 


(102) 


The  inequality  (102)  holds  a  fortiori  if  |  and  q  are  replaced 
by  their  respective  upper  bounds  a  and  b.  Thus,  there  is  finally 
obtained  for  all  points  in  the  given  rectangle,  the  inequality 


Tj/1)  _  Tj/2) 


(104) 


This  series  is  easily  shown,  by  elementary  tests,  to  converge.  Further¬ 
more,  it  is  clear  that  the  quantity  €  can  be  factored  out  of  the 

expression  [”(€  +  m)^  —  nf*'],  so  that  the  series  may  be  written  as  the 
product  of  e  by  a  new  series,  which  is  again  easily  shown  to  converge 
for  al 1  values  of  e.  Designating  this  new  series  by  M,  the  conclusion 
of  the  theorem  is  obtained. 

In  less  precise  language  than  that  used  in  the  formulation  of  the 

theorem,  the  result  obtained  may  be  expressed  by  saying  that,  if  the 
coefficients  f  M  of  the  two  differential  equations  (equation  (92)) 
approximate  each  other  very  closely  in  a  certain  region  and  if  the 
respective  solutions  satisfy  the  same  side  conditions  (equations  (94) 
fl-nri  (95)),  then  the  solutions  approximate  each  other  very  closely 
throughout  the  given  region. 


u 
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In  particular  it  follows  from  this  theorem  that,  if  the  sequence  of 
coefficients  of  the  given  hyperbolic  equation  converges  uniformly,  the 
corresponding  sequence  of  solutions  also  converges  uniformly  to  the 
solution  of  the  limiting  equation. 
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TABLE  I.-  F1  AS  A  FUNCTION  OF  M,  3,  AND  2A 


M 

B 

2A 

F1 

1.00 

0.0000 

0.0000 

00 

1.04 

.2857 

.0120 

1.08 

.4079 

.O338 

1.12 

.5044 

.0606 

1.16 

.5879 

.0910 

15.0172 

1.20 

.6633 

.1242  : 

7.4967 

1.24 

•  7332 

.1595 

4.1109 

1.28 

.7990 

.1964 

2.3598 

1.32 

.8616 

.2346 

1.3609 

1.36 

.9217 

.2738 

.7475 

1.40 

.9798 

.3137 

.3483 

1.44 

I.O361 

.3542 

.0758 

1.48 

1.0911 

.3950 

-.1179 

1.52 

1.1447 

.4362 

-.2606 

1.56 

1.1973 

.4774 

-.3692 

1.60 

1.2490 

.5187 

-.4544 

1.64 

1.2998 

.5600 

-.5231 

1.68 

1.3500 

.6012 

-.5799 

1.72 

1.3994 

.6422 

-.6281 

1.76 

1.4483 

.6829 

-.6699 

1.80 

1.4967 

.7234 

-.7070 

1.84 

1.5445 

.7637 

-.7405 

1.88 

1.5920 

.8035 

-.7713 

1.92 

1.6390 

.8430 

-.8000 

1.96 

1.6857 

.8821 

-.8273 

2.00 

1.7321 

.9208 

-.8533 
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0 

TABLE  II.-  ¥±t  AND  T^2)  AS  FUNCTIONS  OF  2A,  M,  AND  B 


2A 

M 

B 

*1 

r(i) 

r(2) 

0.00 

1.0000 

0 

00 

— 00 

— 00 

.04 

1.0899 

.4334 

84.078 

-2.6449 

-0.1211 

.08 

i.i46o 

.5597 

19.833 

-.9974 

-.0561 

.12 

1.1951 

.6544 

8.1146 

-.4859 

-.0279 

.16 

1.2406 

.7341 

4.0787 

-.2537 

-.0136 

.20 

1.2838 

.8051 

2.2408 

-.1315 

-.0061 

.24 

1.3255 

.8701 

1.2579 

-.0634 

-.0024 

.28 

1.3663 

.9310 

.6737 

-.0257 

-.0007 

.32 

1.4063 

.9887 

.2989 

-.0067 

-.0001 

.36 

1.4457 

1.0441 

.0441 

-.0002 

.0000 

.4o 

1.4848 

1.0976 

-.1375 

-.0022 

-.0000 

.45 

1.5334 

1.1625 

-.3003 

-.0134 

-.0004 

.50 

1.5819 

1.2257 

-.4182 

-.0315 

-.0015 

.55 

1.6303 

1.2876 

-.5077 

-.0548 

-.0036 

.60 

1.6789 

1.3486 

-.5784 

-.0820 

-.0070 

.65 

1.7277 

1.4089 

-.6365 

-.1124 

-.0119 

.70 

1.7768 

1.4687 

—  .6860 

-.1455 

-.0183 

.75 

1.8264 

1.5283 

-.7294 

-.1809 

-.0265 

.80 

1.8765 

1.5878 

-.7686 

-.2184 

-.0364 

.85 

1.9271 

1.6474 

-.8050 

-.2577 

-.0483 

.90 

1.9784 

1.7071 

-.8394 

-.2988 

-.0622 
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2.-  The  strip  <A<~,  o<  a  <  [(1/h)  -  l]  and  the 
le  of  convergence  of  the  development  of  F ^  around  the  point 
°>  A  =  o'71" [(1/h)  -  l]  . 
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0  .  663  .922  1.14  1.35  1.54  ,  1.73 


Figure  4.-  F^  as  a  function  of  M,  B,  and  2a 


